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Abstract 

In the tight-binding approximation we consider multi-channel transmission through a billiard 
coupled to leads. Following Dittes we derive the coupling matrix, the scattering matrix and the 
effective Hamiltonian, but take into account the energy restriction of the conductance band. The 
complex eigenvalues of the effective Hamiltonian define the poles of the scattering matrix. For 
some simple cases, we present exact values for the poles. We derive also the condition for the 
appearance of double poles. 
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I. INTRODUCTION 



In recent years, the ballistic transport through quantum systems has been studied as 
a scattering problem on billiards (microwave cavities) with infinitely high potential walls 
(hard wall approximation). The scattering properties of such billiards are closely related 
to the spectral properties of the corresponding closed billiards The opening of the 

billiards is realized by attaching at least one lead to them. However, for the study of 
the transmission through the billiard two leads are necessary. The fundamental object 
that characterizes the process of quantum scattering is the unitary S-matrix relating the 
amplitudes of incoming waves to the amplitudes of outgoing waves. Provided that the 
properties of the Hamiltonian Hb for the closed billiard are known, one can consider its 
open counterpart and work out the S-matrix formalism by standard methods of the theory of 
quantum scattering . As a result, the S-matrix is expressed in terms of both 

the Hamiltonian Hb and the matrix elements describing the coupling of the billiard states to 
the lead states. The explicit expressions for the coupling matrix elements were formulated 
at first in 1996 by Seba et al. 1G, llll for the case of point contacts of the leads with the 

Q 

billiard. Later, Fyodorov and Sommers p| developed their theory for the connection of the 



billiard with one lead of finite width by using Neumann boundary conditions (see also |2|). 
Recently Dittes |8( considered the same type of an open system and derived the expressions 
for the coupling matrix with Neumann or Dirichlet boundary conditions both by using the 
Green function technique. 

In the present paper, we consider a d-dimensional billiard connected to leads by us- 
ing another approach that is based on the tight-binding model. The motivation for this 
consideration is the following. First of all, the increasing development of fabrication tech- 
niques requires the possibility to perform reliable numerous experiments on ballistic trans- 
port through devices of atomic size 121 1 131 and through molecular devices consisting of very 
few atoms. Secondly, Pichugin et al jljj who applied the formula for the coupling matrix 
derived by Dittes [a], found that their results do not coincide with those of a direct numerical 
computation of the S-matrix poles, above all for the Dirichlet boundary conditions. As we 
will show in the present paper, one of the reasons for this disagreement is that the formal 



continuum approach used by Dittes 



is unbounded in energy and gives zero radiation 



shifts. In electron transmission through electron wires, however, the energy of the electrons 
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is bounded in energy, at least from below. This fact gives rise to radiation shifts of poles 
of the S-matrix which can not be neglected in calculations for concrete systems. Thirdly, 
it is desirable to receive numerical results from an S-matrix computation within the tight- 
binding model in order to compare them with the results of numerical computation of the 
transmission through billiards. To this aim we derive, in the present paper, the coupling 
matrix, the effective Hamiltonian and the poles of the S-matrix within the tight-binding 
model and present some typical numerical results. 

Here, the following remark should be added. The computer simulations solve the 
Schrodinger equation using finite-difference Hamiltonians, i.e. the tight-binding approxima- 
tion. After matching the incoming and outgoing waves with the solutions of the Schrodinger 
equation by implying the boundary conditions at the transverse sections of the leads, the 
conductance of the billiard as well as the scattering wave function can be computed. Today, 
the calculations can be performed with a very high accuracy by using large grids and the 
technique of sparse matrices. The current S-matrix theory is adequate to these computer 
simulations. It is, however, numerically more time consuming because the effective Hamil- 
tonian is not a sparse matrix. Nonetheless, calculations with the effective Hamiltonian are 
useful since they provide another view to the results. In this formalism, the resonant peaks 
of the conductance are related to the poles of the S-matrix that correspond to the eigenval- 
ues of the effective Hamiltonian. It is possible therefore to draw some conclusions on the 
origin of the resonant peaks and on their possible control by means of external parameters. 

In the present paper we will follow, as closely as possible, the Dittes review jsj], even 
in the notations. In the case of microwave or quantum semiconductor billiards, the waves 
are incident to the billiard through (infinitely) long straight waveguides (leads) of a certain 
width. The different channels correspond therefore to different transverse modes of the 
wave propagation within the leads |2|. At a given frequency E (the Fermi energy), we 
enumerate the propagating modes by p = 1,...,M. Thereby, we associate with the lead 
region a continuous set of states \C,p,E > where C specifies the lead number (terminal). In 
the present paper, we consider mostly two leads, the right lead with incident and reflected 
waves and the left lead with outgoing waves as shown in Fig. ^ Our approach can, however, 
be easily generalized to a larger number of leads as done in Sections 4 and 5. 
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II. ONE-DIMENSIONAL TIGHT-BINDING MODEL OF RESONANT TUNNEL- 
ING 



A numerical scheme for the computation of quantum transport through billiards with 
attached leads is mainly based on the finite-difference Schrodinger equation. Implying the 
Ando procedure for the boundary conditions has enabled us to find the transmission 
properties of the billiards from the scattering wave function for any geometry of billiards and 
straight leads. In order to compare the results from such a computation with the S-matrix 
theory, we consider systems projected on a lattice with finite grid. 

As a first example, we consider a simple one-dimensional model for quantum scattering 
and transport. This model is formulated as the tight-binding model (the Anderson model) 

H = -Y,tj\j><j + n+c.c. (i) 

where tj — v l if j = 0, tj — Vr if j — N , and tj = 1 otherwise. This tight-binding model 
presents the simplest case of a one-dimensional box with N sites coupled with left and 
right semi infinite leads via the corresponding coupling constants vl,vr. For a schematic 
representation see Fig. ^ 

The hopping matrix elements tj in the Hamiltonian are proportional to overlapping 
integrals of electron wave functions between adjacent atoms. The model being similar to the 



3Q, 



describes resonant tunneling. 



one-dimensional model with a double barrier structure 
In this case, the values of the coupling coefficients vl,vr play the role of the heights of the 
double barrier structure provided that vl < 1,Vr < 1. The present model gives, however, 
also the possibility to consider the case of strong coupling, vl > 1, Vr > 1. At the left of the 
box we present the solution of the Schrodinger equation 

H\tfj >= E\tfj > (2) 

as 

^ = e ikj + re~ ikj ,j < 1 (3) 
where r is the reflection coefficient with energy 

E(k) = —2 cos k, — 7r < k <ti. (4) 
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FIG. 1: The one-dimensional tight-binding model. The wave lines couple the left and right leads 
with the box containing N points and, correspondingly, N resonant states. 

As will be seen later, the energy E(k) forms the conduction band —2 < E < 2 with finite 
width. At the right of the box we write 

fa = te ikj ,j>N. (5) 

At last, the solution of (J2J) inside the box is 

fa = ae ikj + be~ lk3 ,j = 1, 2, . . . , N. (6) 

Substituting these functions into (J2J) we obtain the following linear equations 

r{e lk -E)+ v L e ik a + v L be~ ik = e ik 
v L r + e ik {e ik - E)a + e~ ik (e~ ik - E)b = -v L 
e ikN (e~ ik - E)a + e~ ikN (e ik - E)b + v R e ik{N+l) t = 

v R e ikN a + v R e~ ikN b + e ik ^ N+1) (e ik -E)t = 0. (7) 

The coefficients a and b can be expressed via the transmission coefficient t as follows 




FIG. 2: The transmission probability versus the wave number of the incident quantum particle. 
The positions of the eigenvalues of the closed billiard (vl = vr= 0) are shown by stars. One can 
see the radiation shifts caused by the coupling of the Id box with the leads. 



1 \„2ikN 



1 a -2ik 



Finally, one obtains 



-2ik^ 



VR 



vr 



for the solution of the system of equations (J7J) where 
.4 



{v L ){vr )e 

v L v r 



2ikN 



e~ 2ik {v L e 2ik 



t = Asin 2 k/A 



-)(v R e 2ik - -). 



(8) 



(9) 



vl v r 

For the particular case vl = vr = 1, we obtain from (j^J) that t = 1 and r = 0. A typical 
resonant transmission through a one- dimensional box is shown in Fig. |21 

The tight-binding model (Jl} demonstrates a few remarkable features. The first one is the 
symmetry of the resonant transmission relative to Vl,r — > 1/vl,r- This symmetry means 
the following: for small as well as for large coupling coefficients, the effective coupling of the 
box to the leads is small. Such a feature was firstly observed in reactions on atomic nuclei, 
see the review J(| , and analytically derived by Dittes et al for an N- level system coupled to 



one open channel With increasing coupling strength, the widths of N — 1 resonance 

states decrease as 1/v in the single-channel case while only one resonance state accumulates 
almost the total sum of the widths. In our case, the iV-level system is coupled to two open 
channels. Correspondingly, with increasing coupling coefficients Vl, Vr the widths of N — 2 
resonance states decrease while the widths of two resonance states increase. We will return 
to this feature below when the poles of the scattering matrix will be considered. The second 
feature is that the heights of the resonant peaks are equal to one only when vl = vr, similar 
o the double barrier resonant structure. This fact was firstly established by Ricco and Azbel 
17| . The radiation shifts of the positions of the resonant peaks relative to the eigenvalues 
of the box 

E n = -2cosk n , k n = 7m/(N + 1), n = 1, 2, . . . , N, (10) 

are the third peculiarity of the tight-binding model. The positions of the eigenvalues (fTUjl 
are shown in Fig. |2]by stars. The shifts and widths of the resonant peaks are symmetrical 
relative to k — > —k, and E — > —E. The last symmetry follows from the invariance of the 
solution of the tight-binding model relative to tj — > —tj. 

III. THE S-MATRIX FOR THE ID TIGHT-BINDING MODEL 

The simplicity of the model allows us to establish the explicit correspondence between 
re analytical results for the transmission amplitudes (jHJ) and the S-matrix approach 



the ana 

ODD 



This Id model was also used in [19( to investigate the width distribution. In 
this approach the scattering system is decomposed into a closed subsystem described by the 
internal Hamiltonian Hb with discrete bound states \ip n >,n = 1,2, N and the continuum 
of external scattering states \E,L > and \E,R> corresponding to the semi infinite left and 
right leads. The Hamiltonian of the two uncoupled subsystems is 

Ho = Hb + Hi + H R , 
H B = 53 E n\ n >< n \, 

n 

H L = J 2 dEE\E,L >< E,L\, H R = J* dEE\E,R>< E,R\], (11) 



7 



where E n are the energies of the bound states of the Id closed billiard, and E denotes the 
energy of the leads. We use the following normalization conditions 



< n \m >= 5 nm 

< E, L\E', L >=< E, R\E', R >= 6(E - E'). (12) 

The couplings between the internal and external subsystems can be incorporated by the 
coupling operator 

V = Y, E f dEV n (E,C)\E,C ><n\+H.C. (13) 

As shown in Fig. ^the closed Id billiard consists of iV sites with energies given by Eq. 
(jlUJl and the corresponding eigenf unctions 



Mj) = J^-rsin(^-),j = l,2,...,N. (11) 



N+l \N + 1. 

These eigenfunctions satisfy the Dirichlet boundary conditions ^n(O) = if) n (N+l) = 0. Since 
the leads are semi-infinite wires the wave functions of the left and right lead are respectively 



1>e,l(J) = \ L | 1 n sinfc(l -j), i/>e,r(J) = \ o l • u\ smfc 0" ~ N )- ( 15 ) 
y 2n\ sin k\ y 2n\ sin k\ 

The energy in the leads corresponding to a single conductance energy band, is defined by 
E(k) = — 2cos/c, —7i < k < 7r. It is easy to see that in the continual limit k — > the 
functions (fT5j) take the form 



17] 



^e,l{x) = J 77-777 sin fcx (16) 
V Ln\k\ 



given in |7[. From (J13)) we have 



V n (E, C) =< E, C\V\n >=< E, C\ £ \j >< j\V\J2 \f >< J> >, C = L, R. 

j j' 

Since as shown in Fig. Q]the coupling matrix elements < J [ V'jj' > are not equaled to zero 
if only j = 0, 1, j' = 1, or j = N, N + 1, j' = N + 1, N, we obtain finally, using (|14|) and 
(JTSJ), the coupling coefficients as 



sink Tin 
V n {E,L) = v L ip E , L (tt)ip n (l) = v L J sin 



vr(iV + l) N + l' 
V n {E, R) = v r MN)1>e,r{N + 1) = W .t nfc ^ sin ^7~TT- (17) 



7r(iV + l) iV+1' 



s 



For the total Hamiltonian 

H = H + V 
the stationary Schrodinger equation reads 

H\ip(E) >= E(k)\ip(E) > . 



(18) 



(19) 



For E(k) different from the eigenvalues E n of the box, the operator (E + zO — Hq)V is well 
defined and equation (flUj) is equivalent to the Lippmann-Schwinger equation 



IV >= |Vo > HE + zO - Hv)- l v\4> > 



(20) 



if the boundary condition of outgoing waves is adopted and 



(£-#o)|Vo >= 0. 
The Lippmann-Schwinger equation (j2Dj) also reads 

|V>= [F(E + iO)]- 1 \i/> >, 



(21) 



(22) 



where 



and 



F(E + iO) = l-(E + iQ- HqY^V 



IVo >= 



\E,L> 




Following to 



(23) 



(24) 



\\E,R> J 

we introduce three projection operators: for the left and right leads 



P c = I dE\E,C>< E,C\ 



and for the billiard 

Pb — \n >< n\ 

n 

with the help of which we can write the scattering wave function (j22j) as 



(25) 



(26) 



|V >= 



Pl|V> 



\ 



Pb\^> 
{Pr\*P> J 



\\i>R> J 



(27) 



Then the coupling operator (fT3]l reads 



PlVP l P l VP b P l VP r 



\ 



( 



V = PbVPl PbVPb PbVPr 
\ P R VP L P r VP b PrVPr ) 

where by using (JT7j) we obtain 



V LB 



V BL V BR 
V Vrb o J 



Vbl = v l Y,M±)^ J dE [l-{E/2f] 1 ' A \n><E 1 L\ 
V BR = vrEMN)^ JdE[l- (£/2) 2 ] 1/4 | 



\n >< E, R\ 



and Vrc = Vq b . Substituting (^Hj) into (J23|) we have 



1 

1 V, 



E^ V LB 



E-H B 




BL 



E^ V BR 



E-H R 



RB 



Using the identity 



-i 



V 



1 -A 
-B 1 -C 
—D 1 



1 1 + ATB AT ATC \ 
TB T TC 
DTB DT 1 + DTC 



one obtains for the inverse matrix F 
( 



F- 1 



where 



and 



1 + E^W£ V LBiE^ V BL e^V LB ^ ^j^Vlb^-^j^Vj 



BR 



1 1 



-V BL 
i i 



D E-H B 

V e-^Vrb^z^Vbl 



l_ 

D 



1 1 



-Vn 



\ 



D E-H B " BR 
^Vrb^ 1 + -^j^Vrb^^j^Vbr J 



1 



1 - BA - CD 



D = 1 — 



E-H 



E v> 



BC 



B C=L,R 



E-H. 



CB- 



C 



From Eqs (J22j) and (J27j) it follows that the wave function in the interior of the billiard 

\^ B >=Q- 1 E v bc \e,c> 

C=L,R 
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where 



Q = E + -H B - £ V BC f+ 1 tfc B . (35) 



Here we used the indentity 



1 :A = A- 1 



l-AB l-BA 

If we substitute the coupling constants (fTTj) into formula (|35|l . it follows that the matrix 
elements of the operator (f33|) can be presented as matrix elements of the effective Hamilto- 



man 



< m\Q\n >= E + 5 mn — < m\H e ff\n > (36) 

where 



1 f 2 Jl - [E x j2Y 

< m\H e ff\n >= E m 8 mn + —V mn / dE x v — = E m S mn - V mn e l . (37) 

2n J -2 E + tO — Ei 



The last expression was obtained by using the formula 

1 1 

iir5(x) + P— 



x + iO x 
where P denotes the principal value integral and 

Vmn = v 2 L iP m (l)iP n (l) + v 2 R iP m (N)ip n (N). (38) 

As can be seen from (|37j) . the lattice approach gives rise to a finite shift of the resonant 
energies 

F m n{E) = -^V mn E (39) 

which is of the same order of magnitude as the width of the resonant peak of the transmission 

lmn(E) = ^V mn /l - (E/2y. (40) 
This energy shift is the main difference between the present tight-binding (lattice) approach 



and the continuum approach by Dittes 



where the shifts are equal to zero. The reason for 



this difference is that the energy is restricted to the conductance band E = —2 cos k. 
TheS- matrixis fla 

Sec = Sea -2m< E^VcbQ^Vbc^C >= r ^ | (41) 
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where r and r' are the reflection coefficients from left to left and from right to right, respec- 
tively, and t and t' are the transmission coefficients from left to right and from right to left. 
Using the definition of the coupling operator (|29|) we can write the transmission coefficient 
of the S-matrix (J41j) as follows 

t = -2iriJ2V m {E,L) <m\Q- 1 \n>V*{E,R). (42) 

mn 

A concrete calculation of the transmission coefficient needs the procedure of inversion of the 
matrix (}36j) . It is therefore more convenient to use a representation by means of the set 
of eigenstates of the effective Hamiltonian (J3~7|) [f], EJ. Using the biorthogonal basis of the 
effective Hamiltonian 

H eff \X) = z x \X), (A|Y) = 5 Xi y, |A) = |A >, (A] =< A|*, (43) 

and the projection operator 

^// = I»(A| (44) 

A 

we obtain 

A E-z x 

Eq. (f4"5|) shows immediately that the eigenvalues of the effective Hamiltonian z\ define 
the poles of the scattering matrix. We underline that the coupling coefficients in the pole 
representation (|43|) have to be calculated by means of the eigenstates |A) of the effective 
Hamiltonian but not with the eigenstates \b > of the closed billiard. The importance of this 
difference is presented in j^j]. 

Let us consider for illustration the limiting case N — 1, the one-sided Anderson model 
[lflj ]. which corresponds to the Id box with a single eigenstate. For simplicity we take 
vl = vr = v. Then the formula for the transmission coefficient © reduces to 

iv 2 sink , , 

1 = 1 Tlk- 46 

cos k — v z e tl: 

On the other hand, from (J37)) and (J38|) we have the effective Hamiltonian as the c-number 
H e ff = z\ = E\ — 2v 2 e lk where for the one-sided dot E\ = 0. Moreover Vi(E,L) = 



Vi(E,R) = V\J^^. Substituting these formulas into (j4^|) we obtain the same formula as 
()46|) with account that E = —2 cos k. An analysis of the S-matrix for the transmission 
through the iV-sided Id box is given in Section 5. 
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FIG. 3: The two-dimensional billiard B attached to three different leads C = 1,2, 3. The coupling 
coefficients x>c between the leads and the billiard are shown by wave lines 

IV. S-MATRIX THEORY FOR THE TRANSMISSION THROUGH BILLIARDS 

Let us consider a d-dimensional billiard specified by the internal eigenstates \b > and 
eigenvalues Eb, 

H B \b >= E b \b > . (47) 

The shape of the billiard is given by the d — 1 surface Q which encloses the internal region of 
the billiard D with the points x G D. The eigenfunctions are < x|6 >= ^b(x). We assume 
that M leads are attached to the billiard. Each lead is a d-dimensional tube with arbitrary 
transverse section uc, C — 1, 2, . . . , M and is semi-infinite along the direction z _L u>c- The 
geometry of the system is illustrated in Fig. 01 for the two-dimensional case and M = 3. It 
allows a separation of variables xi e ujq and z. 

Assuming that the eigenvalues and eigenfunctions for the transverse section of the C-th 
lead are known and denoted by E pc , pc (xj_), we can write the Schrodinger equation for 
the leads in the following manner 

H c \E,C,pc>=E\E,C,pc> ■ (48) 

Here 

E = -2 cos(k Pc ) + E pc , (49) 
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</v( x ) = y 27r|sinfc | shlk ctiz-jc)<l>pc(x±), ( 50 ) 
and jc is the longitudinal position of the attachment of the C-th lead to the billiard. Then 
the Hamiltonian of the uncoupled system consisting of the billiard and the M leads is 

M r 2+E. 



H = Y.E b \b X&I + EE/ dEE[\E,C, Pc >< E } C }Pc \+H.C. (51) 

b C=l Pc j - 2 + e pc 

Similar to (jlHj) let us write the coupling operator as 

^ = EEE / dEV b (E,C,pc)\E,C,p c ><b\ + H.C, (52) 

b C PC J - 2+e pc 

where 

V b (E,C, Pc )=<E,C,pc\V\b> . (53) 

Let be Ac C f2 the areas at which the leads are attached to the billiard. They terminate 
the semi-infinite leads at j z = jc- The shape of Ac is uc, the transverse section of the 
lead. The C-th lead is connected to the billiard through the hopping matrix elements vc, 
as shown in Fig. El by wave lines. Substituting (|5U|) into (p3J), we obtain for the coupling 
matrix elements 



v b (E,c,pc) = Y,^Pc(x)<x\v\y>My) = vJ 1 -^^ E <f> P (y±)A(y±)- (54) 

x,y ' Z7r x x eA c 

It is justified to generalize the one- dimensional case presented in Section 2 to the general 
case with d > 1. In the following, we present some formulas that follow in a straightforward 
manner. Formulas (JHfij) and (J37j) read now 

< b\Q\b' >= E + 5 W - < b\H eff \}f > (55) 



where 



and 



< b\H eff \b' >= E b 5 w -Y,T, W c(b,Pc)W c (b',pc)e l ^c (56) 

C PC 



W c (b,p)=v c E <Mx±)0 P (x ± ). (57) 

The number of channels in each lead is defined by the condition E Pc < E. In the continual 
case the energy is much less than the width of the energy propagation band equaled to 4. 
Therefore we can approximate e lk ~ —E/2 +i. As a result the effective Hamiltonian (|57)|) 
takes the standard form 

H eff = H B -tWW + J (58) 
14 
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FIG. 4: The geometry for the transmission through a rectangular billiard in the tight-binding 
approach. The couplings vq between the leads and the billiard are shown by solid lines 



where Hb is the billiard Hamiltonian the eigenenergies of which are corrected by the radiation 
shifts, and W is the matrix whose elements are given by (JSTj) . The dimension of the matrix 
W is iV x K where N is the number of states in the billiard, K = J2c max(pc'). 

Finally, the S-matrix elements (j4~Tj) are characterized by the channel numbers and read 



< 



C,p c \S\C',p' c , >= 5 CC >5 PP > - 2m < E, C,p c \V C bQ- 1 V B o\E, C',p' c , > . (59) 



For calculation of the S-matrix we can use the set of the eigenstates ipb of the billiard 
or the biorthogonal set of eigenstates ipx of the effective Hamiltonian ()45j) . In the last case 
one can see immediately that the poles of the S-matrix correspond to the eigenvalues of the 
effective Hamiltonian (joT)|) . 



V. SOME APPLICATIONS OF THE GENERAL THEORY 
A. Transmission through a 2d rectangular billiard 

The typical features of the quantum mechanical transmission through a billiard can be 
described by means of a two-dimensional billiard with two attached leads of equal widths 
(Fig. HJ). In the tight-binding formulation x = ao{i,j) where ao is the lattice unit. The 
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eigenf unctions and eigenvalues of the rectangular billiard are 

l/j m ,n(h j) = ^m{i)^n{j) (60) 

E m , n = E m + E n (61) 

where E m and ip m are defined by equations (fTUjl. (fT4*j) for the corresponding numerical size 
of the box N x , N y . Further, the wave functions ()50|) of the leads and their energy ()49|) are 



^E,R, P {hj) = \ L | 1 , I sinfep(z -N x )(f> p (j), (62) 

y ^7T | Sill rCp | 

E = —2 cos A; p + £ p , £ p = -2 cos (^^y) (63) 

where p = 1,2,3,... enumerates the channel number and Nl is the numerical width of 
the leads. We denote the numerical positions of the lead's walls by N\ and N 2 so that 
Nl = N 2 — N — 1. It follows from the geometry of the system that 1 < Ni < N 2 < N y . 
Therefore, the area of intersection between leads and billiard Ac is a straight line of length 
Nl. The eigenfunctions in the transverse sections of the leads have the following form 



Substituting (jfiHj) and ()64j1 into ()54jl we have for the elements of the coupling matrix 
V m JE,L,p) = v L ip m (iy 



sin k p \ 


N 2 

E 


2tt 


sin k p 


N 2 

E 

j=Ni 


2tt 



V m , n (E,R, P ) = VR ^ m {N x )\r-—^ V Mj)M)- (65) 



Here the Latin indexes L, R denote the left and right leads, respectively, as shown in Fig. 
0] The formulas for the effective Hamiltonian (J56j) with (J57j) read 

A 



< m,n\H eff \m ,n >= E m , n 5 mm ,5 nn < - v|Vm(l)Vw(l) £ W L {n,p)W L (n ,p)e*» 

p 

A 

-vl^ m {N x )^ m ,{N x )Y,W R {n,p)W R {n ,p)^, (66) 

p 

with 

W c (n,p)= £ Mj)<PpU)- (67) 

j=Ni 
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The number of channels A is defined by the condition e p < E. 

Let us now consider the correspondence of the formulas obtained to those received in the 
continuum approach by Dittes [8(. First it is necessary to choose, in the last case, some 
characteristic space length. This may be the width of the lead or the size of the billiard. 
Here we choose, as usually in the literature, the former and denote it by d. In the continuum 
approach, the eigenfunctions and the eigenvalues of the rectangular billiard (J5U|) take the 
following form 

4>m,n(x,y) = i)m{x)^n{y)i ^m{x) = \ -SW.(mx/a) (68) 



a 

( 2 



2 2 

m n 



E "" = v ipF + w! (69) 

for Diriclet boundary conditions, where a and b characterize the size of the billiard and 
Eo = h 2 /2md 2 . The eigenfunctions and eigenenergies of the leads are 



^e,l, p {x, y) = J ^j^j s in kx sin npy/ d (70) 

E = E [{kdf + (irp) 2 ] . (71) 
With x = aoi, y = aoj we find the following relations 

1 + N X = a/ao, l + N y = b/ao, l + N L = d/a (72) 

from the comparison of (J60)) and (|62j) with (J68j) and ([70)1 . In the discrete case ip m (j = 0) = 0. 
Therefore, it holds approximately 

, 1pm{l) -^m(O) , v 

for the continuum case, and the coupling matrix elements (}53|) are 

V m , n (E,L,p) = V y' m ^ p (0) (73) 

where 



rV2 

m,n 

and yi, yi are the positions of the lead walls along the y-axis, so that d = y 2 — y\. The same 
expressions were derived in 



It might seem that, for the continuum limit ao — > 0, the radiation shifts go to zero and 
we can use the Weidenmiiller-Dittes approach directly. However as can be seen from (fTTj) . 
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FIG. 5: Two-sided dot coupled to two leads. 

the energy is bounded from below. As a consequence, the principal value integral in the 
matrix elements of the effective Hamiltonian does not vanish. This is in difference to the 
assumption in [s| that E has no limits. 



B. Transmission through a two-sided quantum dot 

A box consisting of a single atom (site) coupled to a left lead and a right one (L and R 
continuous) is the most simple case that gives rise to the Breit-Wigner type formula for the 
transmission amplitude (|46|) shown in Fig. |2] (a). Let us now consider a two sites box that 
gives rise to a 2 x 2 effective Hamiltonian. The properties of such Hamiltonians are studied 
in literature [9j by focusing onto double poles of the S-matrix (branch points in the complex 
energy plane) without relation to a realistic system as well as in relation to laser induced 
structures in atoms . In our present study, we have the possibility to specify the effective 
2x2 Hamiltonian for another specific system to study its properties and to compare the 
results with those of the general study. 

There are different ways to connect the two cite box with the leads as shown in Fig. |SJ 
The cases (a) and (b) are identical, but differ from (c). 

For the case (a) we obtain from (|66J) 

ft//=-r 1+A "i. 



where 



fi -1 + A 



\ = l(v 2 L + v 2 R )e ik , »=\{v\-v\y k . (75) 



The eigenvalues are: 

Zl , 2 = -A ± sjl + p?. (76) 

They define the poles of the S-matrix as shown in Sections 3 and 4. Since the effective 
Hamiltonian is symmetric but not Hermitian we use the biorthogonal basis 0,13 normalized 
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by the condition (m\n) = 5 mn , m = 1, 2, n = 1,2 where {m\ =< m\*. The right eigenstates 
of dHJ) are: 



II >= 



ai > 1 ' -» > '2 >= ° 2 I . (77) 



With the formulas (|56|) . (|57|) the transmission amplitude takes the form 



_ 2iv L v Ry /l - ggg 

The S-matrix has a double pole when two of the eigenvalues ()7(i)l coincide, i.e. when 

E = 0, = 1. (79) 

The energy behavior of the poles (fTHJ) is shown in Fig. |H1 Such a kind of pole behaviour 
was shown in many works based on the general presentation of the effective Hamiltonian 
as a 2 x 2 matrix (see, for example, review [^). The cases (a) and (b) in fig. El {\fi\ > 1) 
correspond to a free crossing of energy levels in the complex plane, while the cases (c) and 
(d) (dashed curves) correspond to the self-avoided crossing. 
The case (c) in Fig. El gives 




Heff=\ , (80) 

\ A l+XJ 
where A is given by (|7H|) . The poles of the scattering matrix are 



2i, 2 = A ± Vl + A 2 (81) 
and the transmission amplitude is given by 



2»tWyi - {E/2f 

(E- Zl )(E-z 2 ) ■ (82) 

A double pole of the S-matrix can be found at 

E = 0, |A| = 1. (83) 

The comparison of (|H2*|) with (fTS|) shows that the way the leads are connected with the 
box, plays an important role for the conductance. In particular, the connection (a) gives rise 
to a transmission zero only at the edges of the energy band E = ±2 while the connection 
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FIG. 6: The energy behavior of the poles for the transmission through the two-sided dot shown 
in Fig. Eta), (a) and (b): = 1 + 0.03. (c) and (d): = 1 — 0.03 (dashed curves). The solid 
curves in (c) and (d) show the case of a double pole, \fi\ = 1. 

(c) leads to t — at E — 0. The energy behavior of the poles (j8Tf is, however, similar to 
that of the poles 1)76)1 and a double pole of the S-matrix appears at E = in both cases. 
Therefore, the transmission is equal to zero in case (c) at the energy where the S-matrix 
has a double pole while this is not so in the cases (a) and (b). This shows clearly that the 
transmission zero for the case (c) is an interference effect. 

C. Transmission through the N-sided Id box 

As a next application, we consider the Id model with N sites presented in Fig. ^in order 
to understand the reduction of the number of transmission peaks by enlarging the coupling 
coefficients. In Fig. |2l the resonant transmission has N = 5 peaks at small coupling 
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V V 
L L 

FIG. 7: Real and imaginary parts of the five poles of the Id chain shown in Fig. ^ versus the 
coupling coefficients vl, vr, E = 1. The chain consists of five sites. 

coefficients vl,vr but N — 2 peaks for large coupling coefficients. 

Using formulas ()37|) and ()38|) . the eigenvalues z^, k — 1, 2, . . . , N of the effective Hamilto- 
nian can be found numerically. In Fig. UJ the real and imaginary parts of the five eigenvalues 
of the effective Hamiltonian (poles of the S-matrix) are shown versus the coupling constants 
vl, -ur for iV = 5 and E = 1. In Fig. UJa, b), the right coupling coefficient vr is chosen to be 
small. In this case, one of the resonance states is broadened with increasing v l and becomes 
shifted beyond the energy band. The incident energy E is tuned to the second energy level 
E 2 = 1 of the box. As a result, this resonance state is broadened. Fig. 0(c, d) demonstrates 
that two resonance states are broadened when both coupling constants are increased. Also 
in this case, the two resonance states are shifted beyond the energy band. Such a nonuniform 
level broadening in the resonance overlapping regime is studied in many different cases by 
using different approaches, see jjj]. It is called resonance trapping The accompanying 
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FIG. 8: The billiard coupled to left and right Id leads at the points jl and j# (a). In (b) 
JL = 3r = jo- For simplicity a rectangular billiard is shown in the figure (although it may be of 
arbitrary shape). The couplings between the billiard and leads vl and vr are shown by wave lines. 

shift in energy appears only when the principal value integral of the matrix elements of the 
effective Hamiltonian is non-vanishing. This is the case in most realistic systems including 
atomic nuclei and atoms [9( and also those considered in the present paper (Fig. (a, c)). 

Thus, the two resonance states do not vanish at strong coupling strength between box 
and leads as it might be concluded from Fig. El The two resonance states go beyond the 
energy band and can therefore contribute to the transmission only via interference with the 
remaining narrow resonant states. This example clearly demonstrates the advantage of the 
effective Hamiltonian approach to the description of transmission. 

D. Transmission through a 2d billiard connected to Id leads 

Let two Id leads be coupled to a 2d billiard at the points jz, (input lead) where jj? (output 
lead) and j = (j x ,j y ), N z = 1 as shown in Fig. |S](a). Substituting the wave functions of the 
Id leads (fT5|) we obtain from Eq. (JHcTJl 

< b\H eff \b' >= E b 5 w + [vlMh)^(h) + v 2 L MiR)^(iR)]e lk . (84) 
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For }l = }r, the equation that defines the poles of the scattering matrix, can be found 
analytically. From (|84|) it follows 



Ei + ^Vi(jo) - e wViCio)V^Cio) (jo)^ (jo) 

^i(jo)^2(jo) E 2 + wV'IQo) - £ ^Mjo^sCjo) 



^i(jo)^s(jo) 



w^(jo)V^Cio) 



£3+^3 2 (jo)-£ 



55) 



where u; = [v L + vj^e is the effective coupling constant. The particular case of a 4 x 4 ef- 
fective Hamiltonian ()85jl was considered in j^J. This determinant can easily be transformed 
to 



IW'Cjo) 



Xi + 1 1 1 
1 x 2 + 1 1 

1 1 £3 + 1 



where 



Ei, — E 



x b 



As a result, the equation for the poles of the S-matrix reads 

K + 4)e-^(j ) 
V £ - E h 



b Xb . 



(86) 



(87) 



E. A 3d billiard connected to a 3d lead 



Consider a 3d billiard that has an arbitrary shape in the x, y plane but is restricted 
in the z-direction by two parallel planes separated by a distance d. In the tight-binding 
approximation, the height of the 3d billiard can be specified by the number N z being equal 
to 1,2,3, .. .. This billiard allows to separate the variables, and it is characterized by the 
eigenvalues and eigenstates of the 3d box Hamiltonian Hb 

H B \b ± ,n z >=(E b± +E nz )\b ± ,n z >, (88) 

where E b± are the transverse eigenenergies of the billiard and 

E nz = -2 cos (^jP^j n z = l,2,...,N z (89) 
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FIG. 9: 3d billiards connected to different 3d leads, (a) The billiard is coupled to iVj, Id leads at 
each point j where iV& is the number of sites of the billiard, (b) The billiard is coupled to one 3d 
lead with the same transverse section as the billiard in 2d. For simplicity, a rectangular billiard is 
shown although it can be of arbitrary shape. The couplings v between billiard and leads are shown 
by wave lines. 

are the longitudinal eigenenergies of the box in ^-direction. The eigenfunctions in the z- 
direction are 

We consider two different types of the connection between the billiard and the leads that 
both are shown in Fig. El In the first case, (a), the Id leads are coupled to the billiard at 
every point j,j z — 1 of the billiard. In the second case, (b), the billiard is coupled to the 3d 
lead the transverse section of which coincides with the shape of the billiard in 2d. 

For the case (a) formula (|33jl reads 

N b i 

Q = E + -H B -Y: V bc f+ „ V CB , (91) 

£7=1 L+ - H C 

where Nj, is the total number of cites of the billiard in the x, y plane. This number is equal 
to the total number of states \b >. The coupling matrix elements are 

/ Sill Aj 

< b,n z \V BL \C >= < >< Ji|VBo|j2 >< h\C >= W'&CjcW^^, (92) 

jl J2 
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where we assume that the coupling between the billiard and each Id lead (with the eigen- 
function (fTSjl) is equaled to v. Substituting (J9~2~|) into (j9*Tj) we obtain, similar to (j37J), 

< b\H eff \b' >= E b 5 w +Y. / dEx < b\V BC \C > — — < C\V CB \b' >= (E b +v 2 e ik )6 w . 

c=1 J hi + ?u — hj\ 

(93) 

That means, for the case (a) in Fig. El the effective Hamiltonian is diagonal with isolated 
poles z b = E b + v 2 e lk , k = acos(—E/2). 

In the case (b) of Fig. the billiard is coupled to one 3d lead, and the total Hamiltonian 



is 



H = H + V 1 

r2+E b 



H = J2(Eb + E nz )\b,n z >< b,n z \+Y, / dEE[\E,b>< E,b\ + \b,R >< E,b\], 

b,n z b J - 2 + E b J 

1^ = EE f +Eh ' dEV bjnz , bl (E)\E,b'><b,n z \+H.C.{94) 



where the eigenfunctions of the lead are 



<j,j z \E,b>= J 1 sinfc 6 (l- (95) 
y 2tt\ sm Kb| 

Here, j z runs along the lead and j runs over the sites in the transverse section of the lead. 
Since the transverse section eigenfunctions of the 3d lead coincide with those of the 2d 
billiard, similar to (JHHJ), we have 

E = -2 cos k b + E b . (96) 
The coupling matrix elements in (|94|) are 

V b , n „ v (E) =< b,n z \V\E,b' >= EE < Hi >< ><3,3z\V\l,l z >< l,l z \E,b' > . (97) 

j,jz u* 

As can be seen from Fig. |§] 

< b\j,j z ><j\V\l,l„ >= vStfSjijrji,,!. 
Substituting into (jHH), we get 



V b>nz ,b'{E) = v ^-^Mj2< b \i >< W >< N *K >= v^-^^5 bb 4 nz {N z ). (98) 
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For the continual case the last expression has to be substituted by ijj'(z = d) 8]. Therefore 
the matrix elements of the effective Hamiltonian are 

sin kb 



{ n 1 r 2+E b 

< b\H eff \b' >= (E b + E nz )5 nz< + v 2 tfj nz (N z )tlj n ,(N z )— / dE' 



E + iO-E' j 
= i(E b + E nz )5 nz< + v 2 e lk ^ nz (N z )^ n ,(N z )]5 w .(99) 

If N g = 1, we obtain an effective Hamiltonian for the case (b) in Fig. that is diagonal in 
the eigen basis of the billiard. This result is similar to that of the case (a). For N z > 1, the 
effective Hamiltonian is also diagonal, however with blocks N z x N z at each diagonal place. 
For N z = 2, with account of (}89|) and (}90|) . the matrix block takes the following form 



H eff=\ 2 2 • 100 



—e ikb E b 1 r '~ • 



Correspondingly, the poles of the S-matrix are 



z nz)b = E b + ^e ik » ± ^1 + ^e 2 ^. (101) 
Using (}96|) . the condition for the double pole can be written down, 

v 2 = 2, E = E b . (102) 



VI. WAVE FUNCTION IN THE INTERIOR OF THE BILLIARD 

The wave function in the interior of the billiard is given by the expression ()H4j) . Using 
the projection operator for the billiard, Pb = J2b \b >< ^1 where \b > are the eigenstates of 
the (closed) billiard, we can rewrite this expression as follows 

<Mx) = E Y.QwME,c, Pc )M*) = £.faMx). (103) 

C,p c bb' b 

The scattering wave function in the interior of the billiard can therefore be expanded in 
the set of eigenfunctions ^t(x) of the Hamiltonian of the closed billiard. The expansion 
coefficients are 

fb = Y,Y,Qbb>Vb>(E,C, Pc ). (104) 

C, PC b> 

The drawback of this representation consists in the fact that the expansion (jlOHj) includes 
the procedure of inversion of the matrix (jo3j) . Similar to (J43j) we can use the set of eigenfunc- 
tions of the effective Hamiltonian for the expansion of the scattering wave function. Using 
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relations (|43|) and f)44j) we can write ()34|) as follows 



(105) 



A 



where the expansion coefficients are 



y A (£,c, Pc ) 

£+-^A 



(106) 



VII. SUMMARY 

In this paper we derived the coupling matrix between a closed billiard and leads attached. 
The knowledge of the coupling matrix gives the explicit expression for the effective Hamil- 
tonian, the S-matrix and the scattering wave function in the interior of the billiard . The 
non hermitian effective Hamiltonian reflects the spectral properties of the closed billiard. 
The eigenvalues of the effective Hamiltonian however are shifted in energy and are complex 
because of the openness of the billiard. 

The theory presented is based on the tight-binding approach. That allows us to estab- 
lish the exact correspondence between the S-matrix theory and numerical calculation of 
the transmission through the billiard that is based on a finite-difference Hamiltonian. The 
present approach can be easily applied to the continual case. The advantage of the effec- 
tive Hamiltonian consists above all in the possibility to interpret numerical results for the 
transmission (Fig. EJ) by means of the poles of the S-matrix (Fig. EJ). The last are the 
eigenvalues of the effective Hamiltonian. It allows us therefore to systematically control the 
transmission through billiards. We presented a few specific examples for which the effective 
Hamiltonian reduces to a complex two by two matrix. 
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